function Sb=GetSFR(SFR,alpha,qav,binary,MminMINE,MminSFR)
% binary means whether your SFR is for binary or all stars, ONLY 0 1
if(alpha<0)
    fprintf('\alpha should be larger than 0')
    Sb=NaN;
end

if(binary>0)
    % whether single stars are counted in SFR or not 
    Sb=(MminMINE/MminSFR)^(1-alpha)*(alpha-2)*SFR /(1+qav) /MminSFR /(alpha-1) ;
else
    Sb=(MminMINE/MminSFR)^(1-alpha)*(alpha-2)*SFR /(3+qav) /MminSFR /(alpha-1) ;
end

end